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ABSTRACT 


The  goal  of  this  project  is  to  develop  methods  for  calculating  molecular 
activation  energies.  The  project  has  been  highly  successful,  as  it  has 
resulted  in  the  development  of  a  state  of  the  art  quantum  chemistry  package  to 
do  these  calculations.  This  package  contains  a  highly  vectorized  Cl  program 
which  appears  to  be  the  fastest  Cl  program  that  is  currently  in  existence. 

Due  to  its  speed  and  efficiency,  the  Cl  program  has  the  potential  for 
widespread  applications  in  many  areas  involving  molecular  energetics.  The 
importance  of  this  program  stems  from  the  fact  that  in  molecular  energy 
calculations,  the  Cl  method  will  yield  the  best  possible  result  that  can  be 
obtained  for  a  given  basis  set.  In  addition  to  the  Cl  program,  this 
comprehensive  quantum  chemistry  package  also  contains  state  of  the  art  SCF, 
geometry  optimization,  and  vibrational  frequency  programs. 
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I.  INTRODUCTION 

The  goal  of  this  project  is  to  calculate  the  amount  of  shocking  or 
jostling  that  a  solid  explosive  can  withstand  before  detonating.  The  ultimate 
objective  is  to  determine,  by  means  of  the  calculations,  how  to  modify  the 
explosive  so  that  it  still  retains  all  of  its  explosive  power  but  becomes  less 
sensitive  to  shock.  These  less  sensitive  but  equally  powerful  explosives 
would  reduce  the  number  of  accidents  that  occur  among  personnel  involved  in 
handling,  storing,  transporting,  and  using  the  explosives. 

The  detonation  process  is  complicated  and  involves  many  factors, 
including  "hot  spots",  propagation  of  acoustic  or  detonation  waves,  inter¬ 
action  (or  scattering)  between  detonation  waves  and  explosive  molecules,  etc. 
However,  the  key  factor  in  an  explosive's  sensitivity  to  shock  is  the 
activation  energy  of  the  molecules  that  make  up  the  explosive.  When  an 
explosive  receives  a  shock,  the  energy  of  the  molecules  contained  in  the 
explosive  is  increased.  In  the  region  of  the  shock,  if  the  energy  of  a 
sufficient  number  of  molecules  is  raised  above  the  activation  energy,  then  the 
detonation  begins. ^Thus  if  one  can  increase  the  activation  energy,  a  bigger 
shock  will  be  required  for  detonation  and  the  explosive  will  be  safer  to 
handle.  Therefore  the  first  step  in  the  study  of  detonations  is  to  determine 
accurate  values  for  the  activation  energy  of  the  explosive  molecules. 

One  can  determine  the  activation  energy  of  an  explosive  molecule  by 
calculating  the  total  energy  of  the  molecule  as  a  function  of  a  reaction 
coordinate  such  as  a  bond  length  or  bond  angle.  At  the  value  of  the  reaction 
coordinate  where  the  total  energy  goes  through  a  maximum,  dissociation,  or 
detonation,  occurs,  and  the  activation  energy  is  the  difference  between  this 
maximum  total  energy  and  the  total  energy  of  the  molecule  in  its  equilibrium 
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state 


[  In  order  for  these  calculated  activation  energies  to  be  sufficiently 

accurate  to  be  useful  In  experimental  and  developmental  work.  It  Is  essential 

that  electron  correlation  effects  be  included  in  the  calculations,  as  electron 

l 

correlation  plays  a  very  crucial  role  in  molecular  energetics.  We  have 
Included  correlation  effects  by  means  of  the  Cl  method,  which  has  the 
capability  of  yielding  the  best  possible  results  that  can  be  obtained  for  a 
given  basis  set. 

Very  substantial  progress  has  been  achieved  in  the  Cl  method,  and  as  a 
result  this  project  is  considered  to  be  highly  successful.  Specifically,  we 
have  developed  a  state  of  the  art  quantum  chemistry  package  that  contains  what 
appears  to  be  the  fastest  Cl  program  that  is  currently  in  existence.  This 
highly  vectori2ed  Cl  program  has  the  potential  for  widespread  applications  in 
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II.  CRAY  QUANTUM  CHEMISTRY  PACKAGE 

The  configuration  interaction  (Cl)  method  is  probably  the  best  known  and 
most  widely  recognized  procedure  for  including  electron  correlation  effects  in 
molecular  energy  calculations.  If  carried  to  completion,  the  Cl  calculation 
will  yield  the  best  possible  result  that  can  be  obtained  for  a  given  basis 
set.  A  complete  or  full  Cl  calculation  Is  one  which  includes  all  of  the 
possible  excitations  of  electrons  from  occupied  to  empty  states.  In  practice 
It  is  often  not  necessary  to  carry  the  Cl  all  the  way  to  completion,  but 
rather  to  converge  the  Cl  to  whatever  degree  of  accuracy  one  desires  for  the 
molecular  energies.  A  converged  Cl  usually  requires  at  least  all  single  and 
double  excitations  and  sometimes  all  triple  and  quadruple  excitations  as 
well.  Since  the  number  of  double,  triple  and  quadruple  excitations  can  be 
quite  large  for  some  energetic  molecules  of  Interest,  the  Cl  calculation  can 
be  very  time-consuming,  and  therefore  one  needs  very  fast  Cl  programs  and  very 
fast  computers  in  order  to  make  the  calculations  feasible.  We  have  obtained 
Dr.  Henry  F.  Schaefer's  Cl  programs*-*®,  which  are  generally  considered  to  be 
the  best  and  fastest  Cl  programs  in  the  world.  In  addition  to  their  speed. 

Dr.  Schaefer's  programs  have  the  capability  of  handling  very  large  numbers  of 
excitations  while  yielding  exact  results  for  the  Cl  energies.  This  last 
feature  is  particularly  important  since  many  other  Cl  programs  are  severely 
limited  in  the  number  of  excitations  that  can  be  included  in  the  Cl 
calculation. 

We  have  made  very  substantial  progress  in  this  area  by  converting  the 
entire  Schaefer  Cl  package  to  the  CRAY.  Not  only  are  these  Cl  programs 
running  on  the  CRAY,  but  also  they  have  been  vectorized  to  take  full  advantage 
of  the  multi -processing  capability  of  the  CRAY.  As  a  result  of  this 
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vector izat ion,  these  Cl  programs  are  now  running  very  fast  and  efficiently. 
Test  calculations  indicate  that  these  vectorized  Cl  programs,  when  run  on  the 
CRAY,  are  probably  the  fastest  Cl  programs  that  are  currently  in  existence. 

In  addition  to  the  Cl  programs,  our  quantum  chemistry  package  contains  a 
state  of  the  art  SCF  program  that  will  handle  large  molecules  and  large  basis 
sets,  and  will  perform  both  RHF  and  ROHF  calculations  on  these  large  systems 
with  d  functions  included  in  the  basis  set.  The  package  also  contains  a  state 
of  the  art  geometry  optimization  program  that  has  the  same  capacity  for 
handling  large  molecules  as  does  the  SCF  program,  and  that  works  well  for  ROHF 
high  spin  states  as  well  as  singlet  states.  This  program  is  probably  the  best 
geometry  optimizer  that  is  currently  available.  In  addition,  our  quantum 
chemistry  package  has  a  new  state  of  the  art  vibrational  frequency  program 
that  yields  highly  accurate  values  for  molecular  vibrational  spectra. 
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III.  ACTIVATION  ENERGY  OF  FLUORINE  AZIDE 


Our  CRAY  quantum  chemistry  package  was  used  to  perform  extensive  Cl 
calculations  of  the  activation  energy  (barrier  height)  of  fluorine  azide 
(FN^),  which  is  a  highly  energetic  molecule.  We  are  grateful  to  Dr.  Harvey 
Michaels  for  a  useful  discussion  in  which  he  provided  us  with  the  theoretical 
geometry  of  FN3  in  its  singlet  ground  state.  Using  this  information  as  a 
guide,  we  optimized  the  geometry  of  the  FN3  singlet  at  the  SCF  level  by  using 
the  Gaussian  82  program**  and  the  6-31G*  basis  set,  and  were  able  to  exactly 
reproduce  Dr.  Michaels'  results.  This  ground  state  geometry  for  FN3  is  shown 
in  Figure  1.  The  values  of  the  bond  lengths  and  bend  angles  are  given  in  the 
first  line  of  Table  I.  As  shown  by  the  figure,  the  ground  state  FN3  molecule 
is  planar  with  the  three  nitrogens  being  nearly  linear.  Having  verified  that 
Gaussian  82  works  well  in  geometry  optimizations  for  the  FN3  singlet,  we  then 
successively  increased  the  central  bond  length  and  optimized  the  other  two 
bond  lengths  and  the  two  bond  angles  at  each  step  of  the  way.  As  expected, 
the  optimization  calculations  show  that  as  the  central  bond  is  stretched,  the 
N-F  separation  decreases,  approaching  its  value  for  a  free  NF  radical  in  the 
closed  shell  singlet  state.  The  other  N-N  distance  also  approaches  its  value 
for  a  free  N2  molecule  in  the  singlet  ground  state.  An  interesting  feature  of 
the  central  bond  stretching  is  that  the  FN3  molecule  becomes  more  bent  as  the 
central  bond  length  is  inceased,  as  indicated  by  the  decreasing  value  of  the 
N-N-N  angle.  The  stretching  of  the  central  bond  was  continued  until  the  total 
energy  went  through  a  maximum  and  began  to  decrease.  All  of  the  above  geometry 
optimization  calculations  were  done  at  the  SCF  level  using  the  6-3 1G*  basis 
set.  The  FN3  optimized  geometry  and  SCF  total  energy  at  each  value  of  the 
central  bond  length  are  given  in  Table  I.  As  shown  by  the  table,  the  SCF 
calculations  give  a  maximum  in  the  total  energy  curve  at  an  N2-NF  separation 


Figure  1.  FN3  ground  state  geometry. 
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Table  I.  Optimized  geometry  and  SCF  total  energy  of  the  FN-j  singlet 
for  fixed  values  of  the  central  bond  length  Rl.  The  top 
line  gives  the  optimized  equilibrium  geometry  of  the  FN-j 
singlet  ground  state.  The  geometry  optimizations  were 
done  at  the  SCF  level  using  the  6-31G*  basis  set.  The 
bond  lengths  and  angles  are  defined  in  Figure  1.  The  SCF 
total  energies,  which  are  given  In  the  last  column,  are 
in  Hartrees,  and  the  bond  lengths  are  in  Angstroms. 
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N2  -108.9439495 

NF (a)  -153.7029725 
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Table  II.  SCF  total  energy  of  the  singlet  ground  state  of  N2  and 
the  open  and  closed  shell  singlet  excited  states  of  NF, 
labeled  NF(a)  and  NF(b)  respectively.  The  SCF 
calculations  were  done  using  the  6-316*  basis  set. 
Energies  are  in  Hartrees. 
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singlet  for  values  of  the  central  bond  length  Rl.  The 
C ISO  and  CISD+OC  calculations  were  done  using  the 
geometries  that  were  optimized  at  the  SCF  level.  Rl  is 
in  Angstroms  and  the  energies  are  in  Hartrees. 


of  approximately  1.59  A,  and  a  barrier  height  of  .47  eV. 

We  also  performed  SCF  calculations  on  the  FN3  decomposition  products, 
which  include  the  N2  singlet  ground  state  and  the  two  NF  singlet  excited 
states.  The  lower  of  these  NF  excited  states,  which  is  an  open  shell  singlet, 
will  be  labeled  NF(a)  while  the  upper  one,  which  is  a  closed  shell  singlet, 
will  be  labeled  NF(b).  Table  II  gives  the  SCF  total  energies  of  these 
states.  According  to  the  SCF  results,  the  separated  products  NF(b)+N2  lie  .41 
eV  below  the  FN3  singlet  ground  state,  while  the  NF(b)  and  NF(a)  states  are 
separated  by  .80  eV. 

After  the  SCF  calculations  were  completed,  we  employed  the  CRAY  Cl 
package  to  perform  Cl  calculations  on  FN3,  using  the  geometries  that  were 
optimized  at  the  SCF  level.  The  Cl  calculations  included  all  single  and 
double  excitations  (CISD)  and  employed  the  6-31G*  basis  set.  No  frozen  core 
was  used  in  these  calculations;  all  of  the  FN3  electrons  were  included  in  the 
Cl,  giving  a  total  of  228,826  configurations  for  the  FN3  singlet.  The 
Davidson  correction  (DC)  for  quadruple  excitations  was  also  computed  and  added 
to  the  CISD  energies.  The  results  of  these  calculations  for  FN3  are  shown  in 
Table  III.  According  to  the  table,  the  CISD  calculation  gives  a  total  energy 
maximum  at  1.71  A  and  a  barrier  height  of  .88  eV,  while  the  CISD+DC 
calculation  gives  a  maximum  at  1.75  A  and  a  barrier  height  of  .95  eV.  These 
results  can  be  compared  with  the  corresponding  SCF  values  of  1.59  A  and  .47 
eV.  Thus  the  CISD  and  CISD+DC  correlation  treatments  produce  a  major  change 
in  the  barrier  height,  causing  it  to  increase  to  approximately  twice  its  value 
at  the  SCF  level.  This  shows  that  Cl  is  essential  to  an  accurate  description 
of  the  FN3  decomposition  process. 


IV.  Cl  FORMALISM 
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We  now  give  a  brief  description  of  the  Cl  formalism.  The  Cl  method  ’ 
is  used  to  solve  the  electronic  Schroedinger  equation: 

Ht  =  Ev  (1) 

where  the  Hamiltonian  H,  in  atomic  units,  is  given  by 
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Here  the  upper-case  Latin  subscripts  represent  the  nuclei,  the  lower-case 
Greek  subscripts  represent  the  electrons,  the  Z's  are  atomic  numbers,  and  the 
r's  (rAB*  rAy‘  ryv)  are  distances  between  the  particles.  The  VNN  term  is 
the  total  nuclear  repulsion  energy. 

In  order  to  solve  the  Schroedinger  equation,  we  write  the  wave  function 
v  as  a  linear  combination  of  symmetry  adapted  configuration  functions  «>s: 


*  =  r  c  ♦  (6) 

s  b  5 

where  the  ♦j's  are  linear  combinations  of  Slater  determinants.  In  a  Cl 


calculation,  one  need  include  only  configurations  of  the  same  symmetry  as  the 
particular  electronic  state  being  investigated  (the  reference  state),  since 
Hamiltonian  matrix  elements  between  configurations  of  different  symmetry  are 
zero.  Hence  the  configuration  functions  are  chosen  to  have  the  symmetry  of 
the  reference  state.  Since  the  total  number  of  configurations  of  all 
symmetries  can  be  very  large,  the  symmetry  adaptation  of  the  configuration 
functions  greatly  reduces  the  size  of  the  Hamiltonian  matrix  and  simplifies 
the  calculation. 


The  coefficients  cs  are  chosen  to  minimize  the  energy 


according  to  the  variation  principle.  The  variation  principle  ensures  that 
the  lowest  calculated  energy  is  an  upper  bound  for  the  exact  lowest  eigenvalue 
of  H  (i.e.,  the  exact  ground  state  energy).  It  also  ensures  that  the  lowest 
calculated  energy  for  each  symmetry  is  an  upper  bound  to  the  exact  lowest 
eigenvalue  of  that  particular  symmetry.  Thus  if  one  is  investigating  an 
excited  state  which  is  the  lowest  of  a  particular  symmetry  type,  one  need 
not  include  any  lower  states  in  the  same  calculation. 

Application  of  the  variation  principle  leads  to  the  well-known  eigenvalue 
equation 

l  CHst-ESs^ct  ■  0  <8) 

which  is  solved,  by  standard  techniques,  for  the  eigenvalues  E  and  coefficients 


ct.  In  the  above  equation,  H  ^  and  S$<.  are  the  Hamiltonian  and  overlap  matrix 
elements,  respectively,  and  are  given  by 


st  "  <4s|H|V 
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The  main  computational  step  in  a  Cl  calculation  is  the  evaluation  of  the 


Hamiltonian  matrix  elements  H$t. 


In  order  to  construct  the  configuration  functions  and  evaluate  the  matrix 


elements,  one  chooses  a  set  of  atomic  basis  functions  xp-  One  then  writes  the 


atomic  orbitals  as  a  linear  combination  of  the  basis  functions 


>.  *  i  Y  U  . 

1  p  P  P1 


where  the  coefficients  U  .  are  determined  by  an  SCF  (self-consistent  field)  or 


similar  calculation.  The  evaluation  of  the  Hamiltonian  matrix  elements  is 


greatly  simplified  if  one  chooses  the  atomic  orbitals  <t>^  to  be  orthonormal. 


The  Hamiltonian  matrix  elements  are  then  given  by 


Hst  =  aijhij  +  Jj|a  biJk^iJk£ 


where 


h1j  *  <♦,<?, )l",!»j<V> 


and  the  coefficients  a*.  and  b*.,,.  are  determined  by  a  projective  reduction 

1  J  1  J  K  K 

.  .  .  .  14  ....  .  .  .  ...  . .  , 


calculation  ,  The  integrals  h„  and  9^ j k £ »  which  are  called  orbital 
integrals,  can  be  obtained  from  analogous  basis  set  integrals,  h^  and 


^ijkz’  the  transformation 


h.  .  =  l  U  .U  .h 
U  pq  Pi  qj  Pq 


q.  ..  „  =  i  U  .U  .U  ,U  .q 
yijk*  pqrs  Pi  qj  rk  siapqrs 


where 


hpq  *  <V?l)lhll*q{;i)> 


••  " 


.•  2.v.v.v;avp:w^vv/aV-.v 


08) 


Ws  *  <*p(7])V(?2,l912l*q('t!)*s(?2)>  ' 

The  transformation  of  the  two-electron  integrals,  Eq.  (16),  can  be  quite  time- 

! 

consuming  and  must  therefore  be  handled  efficiently. 

i 

To  summarize,  the  major  computational  steps  in  a  Cl  calculation  are  the 
following : 

1)  Choose  a  basis  set  and  compute  the  basis  set  integrals 

2)  Determine  the  orbitals  by  means  of  an  SCF  calculation 

i 

3)  Transform  the  basis  set  integrals  to  orbital  integrals 

4)  Choose  and  construct  a  set  of  symmetry  adapted  configuration 
functions  for  the  electronic  state  that  is  being  investigated 

5)  Calculate  the  Hamiltonian  matrix  elements  with  respect  to  the 
conf iguration  functions 

6)  Compute  the  lowest  eigenvalue(s)  and  eigenvector(s)  of  the 
Hamiltonian  matrix 

As  mentioned  previously,  a  full  Cl  calculation  (one  in  which  all  the  configura 
tions  of  a  particular  symmetry  are  included)  yields  the  best  possible  total 
energy  that  can  be  computed  for  a  given  set  of  basis  functions  xp • 


V.  USER  MANUAL 


This  manual  describes  the  input  to  the  programs  in  the  CRAY  Quantum 
Chemistry  Package.  These  programs  provide  the  capability  of  performing 
restricted  Hartree-Fock  (RHF)  closed  shell  SCF  calculations,  RHF  open  shell 
(ROHF)  SCF  calculations.  Cl  calculations,  geometry  optimizations  at  the  RHF 
and  ROHF  levels  of  SCF  theory,  and  harmonic  vibrational  frequency  calculations 
at  the  RHF  and  ROHF  levels. 

Programs  in  this  package  are  limited  to  molecules  whose  geometric 
symmetry  belongs  to  an  Abelian  point  group,  i.e.,  to  D2h  and  subgroups. 
However,  it  is  always  possible  to  employ  less  symmetry  than  a  molecule 
actually  posseses.  For  example,  planar  CH3  (Djh  symmetry)  may  be  run  In  C2v 
point  group  symmetry;  pyramidal  CH3  (C3v)  may  be  run  In  C$  symmetry.  In  this 
way,  any  molecule  may  be  studied. 

Most  of  the  programs  contain  no  fixed  dimensions  except  the  size  of 
available  memory,  which  Is  a  machine  dependent  parameter.  By  increasing  this 
dimension,  one  may  apply  these  methods  to  increasingly  larger  systems. 


OVERVIEW 


The  SCF,  Cl,  geometry  optimization,  and  force  constant  evaluation 
procedures  consist  of  several  discreet  programs,  as  contrasted  with  one  huge 
program  module.  Each  of  these  programs  is  summarized  below. 

INPUT 

Program  INPUT  takes  a  simple  formatted  input  deck  and  generates  an 
unformatted  direct  access  file,  FTN010.  This  is  a  master  data  file  containing 
such  information  as  the  basis  set,  the  AG-SO  transformation  matrix,  and  the 
atomic  charges.  INPUT  need  be  run  only  once  for  a  given  system. 

INTS 

This  program  generates  integrals  over  symmetry-adapted  atomic-centered 
Gaussian  basis  functions.  These  include  the  one-electron  (S,  T,  V)  integrals, 
as  well  as  the  two-electron  terms,  <yv|xo>. 

SCF 

The  SCF  program  solves  the  Hartree-Fock-Roothan  equations  for  the  self- 
consistent-field  energy  and  MO  coefficients.  The  MO  coeff icients,  or  SCF 
vector,  are  stored  on  the  master  file,  FTN010. 

DRT 

This  program  generates  the  distinct  row  table  (DRT)  which  is  necessary 
for  the  graphical  unitary  group  approach  (GUGA)  configuration  interaction  (Cl) 
calculation.  Information  contained  In  the  DRT  is  also  required  for  the 
vibrational  frequency  calculation. 


TRANS 


This  program  performs  the  transformation  of  integrals  from  the  atomic- 
orbital  to  the  molecular -orbital  basis  using  the  SCF  vector. 

CISORT 

Program  CISORT  creates  the  3  and  4-external  element  supermatrices  for  the 
GU6A  Cl  program. 

Cl 

This  program  performs  Cl  calculations  using  the  GUGA  method. 

AODER 

Program  AOOER  evaluates  the  first  and  second  derivatives  of  the  AO 
integrals,  with  respect  to  atomic  displacements.  The  derivative  integrals  are 
multiplied  by  the  appropriate  density  matrix  element  to  yield  the  atomic- 
orbital  contribution  to  the  energy  gradients. 

CPHF 

Solution  of  the  coupled-perturbed  Hartree-Fock  equations  is  performed  by 
this  program.  The  first-order  changes  in  the  SCF  coefficients  are  required 
for  an  SCF  force  constant  evaluation. 

OPTIM 

Given  an  energy  gradient  (at  any  level  of  theory)  OPTIM  updates  the 
molecular  geometry  In  the  direction  of  the  energy  minimum.  When  the  energy 
gradient  falls  below  an  Indicated  threshold,  OPTIM  terminates  the  job. 


NORMCO 


Given  a  force  constant  matrix,  this  routine  performs  the  diagonal Ization 
In  mass  weighted  coordinates,  yielding  the  harmonic  vibrational  frequencies 
and  the  normal  modes. 


INPUT  PARAMETERS 


f 

I 

Input  to  each  program  comes  from  both  a  formatted  input  deck  (unit  5)  and  j 

an  unformatted,  direct-access  file  (FTN010),  which  is  generated  by  program 
INPUT.  Each  program  rewinds  unit  5  and  searches  for  a  10  character  flag  which  I 

indicates  the  start  of  data  for  that  particular  program.  This  allows  the 
input  to  the  various  programs  to  be  entered  in  any  order.  Unrecognized  flags 
are  ignored  and  may  be  used  to  "comment"  the  input  deck.  If  input  to  a 
program  is  missing.  It  will  result  in  an  abort,  or  in  the  selection  of  default 
values.  If  possible. 


Program  INPUT 

It  is  the  function  of  INPUT  to  generate  the  unformatted,  direct- 
access  file  FTN010,  which  is  read  by  all  subsequent  programs.  The  input 
to  INPUT  is  designed  to  be  compatible  with  the  HONDO  program  (Dupuis, 
King,  Rys.,  J.  Chem.  Phys.,  65,  111,  1976).  However,  not  all  the  HONDO 
parameters  are  required  here.  Hence,  a  number  of  dummy  parametrs  must 
be  entered,  for  compatibility.  The  following  input  description  is  taken 


from 

the  comment 

cards  in  the  HONDO  program. 

-0- 

FLAG 

ft  INPUT  ftft 

START  OF  INPUT  FOR  INPUT 

-1- 

TITLE 

FORMAT (10A8) 

-2- 

TIMLIM , NPRINT , ITOL , ICUT  FORMAT ( F10 . 0 , 31 3 ) 

OPTIONS  FOR  THE  JOB 

—  TIMLIM  =  INTERNAL  TIME  LIMIT  TO  STOP  THE  JOB. 

NOT  USED  HERE. 

—  NPRINT  =  0  ...  NORMAL  PRINTING 

=  1  ...  EXTRA  PRINTING  FOR  BASIS  SET+SYMMETRY 

=  2  ...  EXTRA  PRINTING  FOR  MO  GUESS  ROUTINES 

=  3  ...  PRINT  OUT  IE- INTEGRALS 

=  4  ...  PRINT  OUT  2E-INTEGRALS 

=  5  ...  PRINT  OUT  SCF  DATA  FOR  EACH  CYCLE 

=  6  ...  LIMITED  PRINT  FOR  TSO  (72  CHAR. /LINE) 

=  7  ...  PUNCH  OUT  BASIS  SET  INPUT  DATA  +  MO'S 

---  ITOL  :  PRODUCTS  OF  PRIMITIVES  WHOSE 

PREEXPONENTIAL  FACTOR  EXP(-G)  IS  LESS 
THAN  10. 0**( -ITOL)  ARE  SKIPPED. 

DEFAULT  VALUE  :  10.0**(-9) 

ICUT  =  INTEGRALS  LESS  THAN  10.0**(-ICUT)  IN 
ABSOLUTE  VALUE  ARE  NOT  STORED  ON  THE 
INTEGRAL  FILE. 

DEFAULT  VALUE  :  10.0** (-9) 


-3-  SCHOENFLIES  SYMBOL  OF  THE  SYMMETRY  POINT  GROUP. 
PROGRAM  IS  LIMITED  TO  D2H  AND  SUBGROUPS. 

FORMAT (A5, 15) 

POSSIBLE  SYMBOLS:  'Cl 

'CS 
'Cl 
'CNH 
'CNV 


I 

I 


I  =0RDER  OF  PRINCIPAL  AXIS 


-4-  ORIENTATION  OF  LOCAL  SYMMETRY  FRAME 

XI,  Yl,  Zl,  X2 ,  Y2,  Z2  F0RMAT(6F10. 5) 

POINT  1  =  (XI , Yl ,Z1 )  POINT  2  =  (X2,Y2,Z2)  TWO  DISTINCT  POINTS 
USE  CARDS  4  AND  5  TO  SPECIFY  NONSTANDARD  ORIENTATION  OF 
SYMMETRY  ELEMENTS. 

DEFAULT  OPTION  IS  TO  LET  SYMMETRY  FRAME  COINCIDE  WITH  MASTER  FRAME. 
TO  ELECT  DEFAULT  OPTION  FOR  ANY  GROUP  OTHER  THAN  Cl  LEAVE  CARD  4 
BLANK,  AND  OMIT  CARD  5. 


GROUP:  SYMMETRY  ELEMENTS  WITH  RESPECT  TO  LOCAL  X,Y,Z  FRAME: 


NO  SYMMETRY.  OMIT  CARDS  4  AND  3. 

REFLECTION  PLANE  =  (X,Y) 

CENTER  OF  INVERSION  IS  ORIGIN,  X=Y=Z=0 
N-FOLD  AXIS  =  Z,  SIGMA-H  PLANE  =  (X,Y) 

N-FOLD  AXIS  =  Z,  SIGMA-V  PLANE  =  (X,Z) 

N-FOLD  AXIS  =  Z,  2-FOLD  AXIS  =  X,  SIGMA-H  PLANE  =  (X,Y) 


MEANING  OF  POINTS  1  AND  2: 

FOR  Cl  GROUP,  NO  CARDS  4  OR  5,  GO  TO  CARD  -6-. 

FOR  CS  GROUP,  ANY  TWO  POINTS  IN  THE  SYMMETRY  PLANE. 

FOR  Cl  GROUP,  POINT  1  =  CENTER  OF  INVERTION.  (IGNORE  POINT  2) 
FOR  ALL  OTHER  GROUPS  POINTS  1  AND  2  MAY  BE  ANY  TWO  POINTS  ON 
LOCAL  Z  AXIS.  (SEE  TABLE  ABOVE) 


-5-  ORIENTATION  OF  SYMMETRY  FRAME  CONTINUED. 

X3,  Y3,  Z3,  DIRECTIONAL  PARAMETER  F0RMAT(3F10.5,A8) 
POINT  3  =  (X3,Y3,Z3)  MUST  NOT  BE  COL1NEAR  WITH  POINTS 

1  AND  2  ON  CARD  4. 

OMIT  CARD  5  IF  DEFAULT  OPTION  IS  ELECTED  ON  CARD  4. 


FOR  Cl  GROUP,  OMIT  CARD  3. 

FOR  CS  GROUP,  POINT  3  =  ANY  POINT  IN  SYM  PLANE.  IGNORE  PARAMETER. 
FOR  Cl  GROUP,  OMIT  CARD  5. 

FOR  OTHER  GROUPS,  INPUT  POINT  3  AND  DIRECTIONAL  PARAMETER. 
POSSIBLE  VALUES  OF  DIRECTIONAL  PARAMETER: 

'PARALLEL'  POINT  3  LIES  ON  LOCAL  X  AXIS.  (DEFAULT  VALUE) 
'NORMAL  '  POINT  3  LIES  ON  LOCAL  Y  AXIS. 


COMMENT:  THE  MOST  COMMON  USE  OF  THE  'NORMAL  '  OPTION  IS  WITH 

DNH  WITH  ODD  N,  E.G.  STAGGERED  ETHANE.  THEN  THE  LOCAL 
PLANE  IS  A  SIGMA-V  PLANE  NORMAL  TO  THE  X  AXIS. 


-6-  CHARGE  ,  MULTIPLICITY  .UNITS.  FORMAT(3I5) 


*******  NONE  OF  THESE  PARAMETERS  IS  USED  HERE  **** 
*******  CHARGE,  MULTIPLICITY  ENTERED  LATER;  **** 
*******  COORDINATES  ARE  ALWAYS  IN  ATOMIC  UNITS  (BOHRS)  **** 


CHARGE  =  +1  FOR  NH4+,  -1  FOR  F-  ION,  FOR  EXAMPLE. 
MULTIPLICITY  =  1  FOR  CLOSED- SHELL  MOLECULE. 

I UNIT  >  0  FOR  ANGSTROMS. 

I UNIT  <  OR  =  0  FOR  ATOMIC  UNITS. 
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-7-  FIRST  UNIQUE  CENTER  FORMAT* 2A4 , 2X, F5. 0, 3F20. 10) 
NAME,ZNUC,X, Y,Z 

NAME  =  ATOMIC  NAME  (USED  ONLY  FOR  PRINT  OUT). 
ZNUC  =  NUCLEAR  CHARGE. 

X, Y, Z  =  COORDINATES  IN  THE  MASTER  FRAME  IN  BOHR. 


-8-  ISHELL,  ITYPE,  IGAUSS 
FORMAT* 15 , IX, A4, 15) 

ISHELL  =  SHELL  NUMBER  ON  THIS  ATOM. 

FOR  EACH  ATOM,  ISHELL  TAKES  VALUES  1,2,... 

ITYPE  =  SHELL  TYPE  ='  S'  FOR  S  SHELLS 

'  P'  FOR  P  SHELLS 
'  D'  FOR  D  SHELLS 

IGAUSS  =  NUMBER  OF  CONTRACTED  PRIMITIVES  IN  SHELL 
***  SPECIAL  OPTION: 

THE  FOLLOWING  CODES  WILL  CAUSE  BASIS  SETS  TO 

BE  READ  DIRECTLY  FROM  THE  BERKELEY  BASIS  SET  LIBRARY. 

USE  THESE  IN  PLACE  OF  ISHELL,  ITYPR ,  IGAUSS. 

GET  DZ  -  USE  HUZINAGA/DUNNING  (9S5P/4S2P)  BASIS  FOR  H-F 

( 1 1S7P/6S4P)  BASIS  AL-CL 
GET  DZP  -  AS  ABOVE  WITH  POLARIZATION. 

ALSO  AUGMENTED  WACHTERS  BASIS  FOR  SE-ZN. 

GET  TZ  -  USE  HUZINGA/DUNNING  (9S5P/5S3P)  BASIS  H-F 

GET  WACHTERS  -  USE  WACHTERS  BASIS  SETS  FOR  TRANSITION  METALS 

GET  ST0-3G  -  USE  ST0-3G  BASIS  (RIPPED  OFF  FROM  HONDO)  H-AR 

IT  IS  POSSIBLE  TO  MIX  "GET"  CARDS  WITH  ORDINARY  BASIS  SETS 
ON  A  GIVEN  CENTER  (SEE  EXAMPLES) 


IF  SPECIAL  "GET"  TYPE  BASIS  FUNCTION,  NO  CARD  9.  GO  BACK  TO  -8- 
-9-  IG,  A,  Cl  FORM AT ( 1 5, El 5. 9, E20. 10) 

IG  =  PRIMITIVE  NUMBER  IN  THE  SHELL  ISHELL 

FOR  EACH  SHELL,  IG  TAKES  VALUES  1 , 2 ,..., IGAUSS. 

A  =  GAUSSIAN  EXPONENTIAL  PARAMETER  OF  PRIMITIVE 
Cl  =  CONTRACTION  COEFFICIENT 

.  THERE  ARE  IGAUSS  CARDS,  ONE  FOR  EACH  PRIMITIVE  OF  SHELL. 

.  THEN  DATA  FOR  THE  NEXT  SHELL,  RETURN  TO  CARD-8-. 


-10-  A  BLANK  CARD  ENDS  LIST  OF  SHELLS  CENTERED  ON  THIS  ATOM. 
.  THEN  NEXT  UNIQUE  CENTER.  GO  BACK  TO  CARD-7- 


-11-  A  BLANK  CARD  ENDS  THE  LIST  OF  THE  UNIQUE  CENTERS. 


Program  INTS 


INTS  receives  all  data  from  file  FTN010;  no  further  input  is 
required.  However,  one  may  optionally  input  new  geometric  coordinates 
if  desired. 


-0-  FLAG  $NEWGE0M  FORMAT(AIO) 

START  OF  NEW  GEOMETRY  FOR  INTS 

if  absent  coordinates  will  be  read  from  FTN010 
if  present  INTS  reads  the  following  cards 


-1-  CHG , X , Y , 2  F0RMAT(4F20. 10) 

charge,  x,  y,  z,  coordinates  for  I-th  atom 
REPEAT  FOR  ALL  ATOMS 

NOTE:  While  INPUT  requires  only  unique  atomic  centers, 
SNEWGEOM  requires  ALL  coordinates 


n 


Program  SCF 


-0-  FLAG 

#  SCF  #### 

FORMAT(AIO) 

START  OF  SCF 

INPUT 

-1-  ALABEL 

FORMAT (80A1) 

-2-  IPRCT,  I SCF,  IOPEN,  INFLG,  IVECT,  IORDER,  PRINT,  IDEBUG,  ISAVE, 
ITRAS,  MAXNO ,  ISTO,  MICMX,  NCOR,  IDFIX 

PARAMETERS  FOR  THE  SCF  CALCULATION  FORMAT(15I4) 

IPRCT  :  NUMBER  OF  ITERATIONS  BEFORE  EXTRAPOLATION. 

IPRCT=0  MEANS  NO  EXTRAPOLATION.  (RECOMMENDED) 

ISCF  :  SCF  IS  CONVERGED  WHEN  RMS  CHANGE  IN  THE  DENSITY  MATRIX 

IS  <  10**( -ISCF) .  DEFAULT=6 ,  BUT  9  RECOMMENDED  FOR 
GRADIENT  CALCULATIONS. 

IOPEN  :  =  0  FOR  CLOSED-SHELL  CALCULATION 

=  1  FOR  OPEN-SHELL  RHF  CALCULATION 

INFLG  :  =  0  TO  GENERATE  NEV  SCF  VECTOR  FROM  BARE-NUCLEUS  GUESS 

:  =  1  TO  READ  OLD  VECTOR  FROM  FTNOIO 

**  OLD  VECTOR  MUST  EXIST  IF  INFLG=1  ** 

IVECT  :  =  0  RE-ORDER  MO'S  BY  EIGENVALUE 

:  =  -4  DO  NOT  RE-ORDER  VECTOR 

IORDER  :  =  0  USE  OLD  VECTOR  IN  PRESENT  ORDER 

:  <>  0  RE-ORDER  MO'S  AS  INDICATED  ON  NEXT  CARDS 

PRINT  :  PRINT  OPTION  FOR  DEBUGGING 

IDEBUG  :  MORE  PRINT  OPTION  FOR  DEBUGGING 

ISAVE:  NOT  USED 
ITRAS:  NOT  USED 

MAXNO  :  MAXIMUM  NUMBER  OF  ITERATIONS.  DEFAULT=40. 

ISTO  :  10**( -ISTO)  IS  THE  THRESHOLD  FOR  DROPPING  INTEGRALS  IN 
THE  P  AND  K  TAPE  CONSTRUCTION.  DEFAULT=20. 

MICMX  :  NUMBER  OF  MICRO-ITERATIONS  (NOT  ENABLED) 

NCOR  :  UNIT  NUMBER  TO  READ  RELATIVISTIC  CORE  INTEGRALS 
(NOT  ENABLED) 

IDFIX  :  =  0  COMPUTE  OPTIMUM  DAMPING  FACTOR  EACH  ITERATION 
:  <>  0  USE  INPUT  DAMPING  FACTOR  FOR  ALL  ITERATIONS 
**  WHEN  CONVERGENCE  IS  <  10**(-5)  DAMPING  IS  ALWAYS 
TURNED  OFF 
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-3-  IF  I ORDER  =  0  NO  CARD  -3-;  PROCEED  TO  -4- 
IF  IORDER  =  1  THEN  NEW  ORDER  OF  THE  MO'S 
FORMAT( 1714) 

STARTING  EACH  SYMMETRY  BLOCK  ON  A  NEW  CARD,  GIVE  THE  NEW 
POSITION  OF  EACH  COLUMN  OF  THE  SCF  VECTOR.  IF  NO  RE-ORDERING 
IS  REQUIRED  FOR  A  GIVEN  SYMMETRY  BLOCK,  BUT  IORDER=l, 

A  UNIT  MATRIX  MUST  BE  INPUT.  (SEE  "REORDERING  THE  SCF  VECTOR"). 


-4-  NCLOS,  NOPEN  FORMAT(2I5) 

ONE  CARD  FOR  EACH  SYMMETRY  BLOCK 

NCLOS  :  NUMBER  OF  CLOSED-SHELL  ORBITALS  IN  THE  SYMMETRY  BLOCK 
NOPEN  :  NUMBER  OF  OPEN- SHELL  OBITALS  IN  THE  SYMMETRY  BLOCK 

**  NOTES:  "MISSING"  SYMMETRY  BLOCKS  (EG.  A,  OF  DZ  H„0)  MUST 
APPEAR  WITH  NCL0S=N0PEN=0 

SYMMETRY  BLOCKS  ARE  ASSUMED  TO  BE  IN  COTTON'S  ORDER, 
F.A.  COTTON,  CHEMICAL  APPLICATIONS  OF  GROUP  THEORY, 
WILEY- INTERSCIENCE,  NEW  YORK,  1971. 


-5-  IF  IOPEN  =  0  NO  CARD  -5-;  PROCEED  TO  -6- 

IF  IOPEN  =  1  THEN  NOP  NUMBER  OF  CARDS,  NOP  =  N*(N+l)/2, 

WHERE  N  =  #  OF  SYMMETRY  TYPES  WITH  OPEN  SHELLS 

ALPHA,  BETA  F0RMAT(2F20. 10) 

COUPLING  COEFFICIENTS  FOR  RHF  OPEN-SHELL  STATES  (SEE  OHIO  STATE 
UNIVERSITY  THEORETICAL  CHEMISTRY  GROUP  REPORT  «  101) 

MOST  COMMON  TYPES  ARE: 

HIGH  SPIN  COUPLED  ELECTRONS  - 

ALPHA  =0.0,  BETA  =  -1.0  FOR  ALL  NOP  CARDS 


OPEN-SHELL  SINGLET  WAVEFUNCTION  - 

(  2  UNPAIRED  ELECTRONS  IN  DIFFERENT  SYMMETRY  BLOCKS) 
ALPHA  =  0.0  BETA  =  -1.0  FOR  CARD  1 

ALPHA  =0.0  BETA  =  3.0  FOR  CARD  2 

ALPHA  =  0.0  BETA  =  -1.0  FOR  CARD  3 


-7-  DAMP,  SHIFTL  FORMAT( 2F20. 10) 

DAMP  :  DAMPING  FOR  ORBITAL  CHANGES  (  0.0  RECOMMENDED) 

SHIFTL  :  LEVEL  SHIFT  FOR  VIRTUAL  ORBITALS  (1.0  RECOMMENDED) 


In  the  case  of  gradient  calculation,  subroutine  TFOCK  of  the  SCF 
program  requires  additional  input. 


-0-  FUG  #  TFOCK  ##  FORMAT(AIO) 

START  OF  TFOCK  INPUT  (FOR  GRADIENTS) 

-1-  SCFTYP,  CALTYP,  DERTYPE  FORMAT( 3A6) 

SCFTYP  -  INDICATES  THE  TYPE  OF  FORMALISM  USED  USED  TO 
EVALUATE  THE  GRADIENT 

"GRSCF  "  SHOULD  BE  USED  WITH  THIS  SET  OF  PROGRAMS, 
WHICH  EMPLOY  THE  GENERAL  SCF  FORMALISM 
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Hucnp  h  rnn  utru  CDTU  Alin 

"CLSCF  H  FOR  CLOSED-SHELLS  MIGHT  BE  REQUIRED  WITH  EARLIER 
VERSIONS  OF  THE  PROGRAMS 

CALTYP  -  INDICATES  THE  TYPE  OF  ENERGY  GRADIENT  TO  BE  PERFORMED 
"SCF  “  FOR  SCF  ENERGY  GRADIENTS 
“MCSCf  "  FOR  MCSCF  GRADIENTS  (NOT  ENABLED) 

“Cl  "  FOR  Cl  GRADIENTS  (NOT  ENABLED) 

DERTYP  -  INDICATES  THE  ORDER  OF  THE  GRADIENT  CALCULATION 
“FIRST  "  FOR  1ST  DERIVATIVES 
“SECOND"  FOR  2ND  DERIVATIVES 


-2-  IPRINT  FORMAT( 15) 

PRINT  OPTION  FOR  DEBUGGING 
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Program  DRT 

The  Distinct  Row  Table  program  generates  a  series  of  arrays  required  by 
the  TRANS,  Cl  and  CP HF  programs. 


-0- 

FLAG 

#  DRT  #### 

FORMAT(AIO) 

START  OF  INPUT 

FOR  DRT 

-1- 

ALABLE 

FORMAT ( 80 Al) 

-2-  OPTIONS  CARD  FOR  Cl  CALCULATION.  FOR  SCF  DERIVATIVES, 

LEAVE  THIS  CARD  BLANK. 

OPTION(l) ,  OPTION (2)  F0RMAT(2I5) 

OPT  ION ( 1)  -  PRINTING  OPTION  (POWERS  OF  TWO  TO  TURN  ON  INDIVIDUAL  ONES) 

1  PRINT  DRT 

2  PRINT  EXTERNAL  WEIGHT  ARRAYS 
4  PRINT  INTEGRAL  POINTER  ARRAYS 

OPTION(2)  -  EXCITATION  LEVEL  FOR  EXCITATIONS  INTO  VIRTUALS  (DEFAULT  2) 


-3-  NSYM,  NBF  FORMAT(2I5) 

NSYM  -  NUMBER  OF  ELEMENTS  IN  THIS  POINT  GROUP. 

THIS  NUMBER  MUST  INCLUDE  “MISSING"  SYMMETRY  TYPES 
(NSYM=4  FOR  DZ  H20,  EVEN  THOUGH  THERE  ARE  NO  A2  ORBITALS) 

NBF  -  TOTAL  NUMBER  OF  BASIS  FUNCTIONS 


-4-  ICOUNT,  KEY,  ISYM  FREE  FORMAT 

ICOUNT  -  REPEAT  COUNT  FOR  NUMBER  OF  EQUIVALENT 
ORBITALS  OF  THIS  TYPE 

KEY  -  ORBITAL  CODE 

FOR  SCF  FORCE  CONSTANTS 

VAl in  TYPF^  ABF • 

DOC  -  DOUBLY  OCCUPIED  ORBITAL 

UOC  -  UNOCCUPIED  ORBITAL 

ALP  -  ALPHA  (HIGH  SPIN)  OCCUPIED  ORBITAL 

SPE  -  SPECIAL  ORBITALS  DEFINED  LATER 

FOR  Cl  ONE  MAY  ALSO  HAVE: 

FZC  -  FROZEN  CORE  ORBITAL 
FZV  -  FROZEN  VIRTUAL  ORBITAL 

ISYM  -  NUMBER  SPECIFYING  THE  SYMMETRY  TYPE  OF  THIS  ORBITAL. 
VALUES  RANGE  FROM  1  TO  NSYM. 

REPEAT  -4-  FOR  ALL  ORBTALS 


26 


-5-  SPECIAL  CODE  FREE  FORMAT 

IF  SPE  APPEARS  AS  AN  ORBITAL  CODE,  THE  TYPE  OF  SPECIAL 
ORBITAL  MUST  BE  IDENTIFIED.  CURRENTLY,  THE  ONLY  VALID 
CODE  IS: 

OSS 

INDICATING  AN  OPEN-SHELL  SINGLET  CALCULATION 
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Program  TRANS  requires  no  Input. 
Program  Cl  SORT  requires  no  input. 
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Program  Cl 


-0-  FLAG  #  Cl  #####  FORMAT(AIO) 


-1-  IGUESS,  MXITER,  IROOTI,  NROOTS  F0RMAT(4I5) 

IGUESS  =  0,  UNIT  VECTOR  BASED  ON  LOWEST  DIAGONAL  ELEMENT 
>  0,  READ  IGUESS  VECTORS  FROM  FILE  12  TO  START. 

MXITER  =  MAXIMUM  ITERATIONS  PER  ROOT  SOUGHT  (10  WITH  2  ROOTS 
GIVES  20  TOTAL). 

IROOTI  =  FIRST  ROOT  SOUGHT,  DEFAULTS  TO  1  WITHOUT  STARTING 
VECTORS  FOR  LOWER  ROOTS. 

NROOTS  =  NUMBER  OF  TOTAL  ROOTS  TO  BE  SOUGHT. 
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Program  AODER 


In  the  absence  of  input,  AODER  will  default  to  the  same  calculation 
as  indicated  for  TFOCK.  Alternatively,  the  type  of  calculation  may  be 
overridden,  and  new  options  specified. 


-0-  FUG  #  NEVDER  #  FORMAT  (A10) 

START  OF  INPUT  FOR  AODER 


-1-  CALTYP ,  DERTYP,  I  PRINT  FORMAT ( 2 A6 , 15) 

CALTYP,  DERTYP  -  SAME  AS  FOR  TFOCK 
1PRINT  -  PRINT  OPTION 


-2-  NOSYM,  INTOUT,  INFOUT  F0RMAT(A5,2I5) 

NOSYM  -  FOR  DEBUGGING,  IF  NOSYM  =  "NOSYM",  DISABLE  SYMMETRY 
IN  THE  DERIVATIVE  EVALUATION. 

-  ELSE  USE  AVAIUBLE  SYMMETRY. 

USE  NOSYM  =  " 

INTOUT  -  <>  0  WRITE  OUT  DERIVATIVE  INTEGRALS 

INFOUT  -  <>  0  WRITE  OUT  TWO-ELECTRON  CONTRIBUTION  TO  DERIVATIVE 
FOCK  MATRIX  (FOR  CPHF): 


Program  CPHF 

The  CPHF  program  will  use  default  values  if  no  input  is  found.  The 
input  parameters  indicated  here  are  optional 


-0-  FLAG  #  CPHF  ###  FORMAT(AIO) 

START  OF  INPUT  FOR  CPHF 


-1-  IPRINT,  ICNVG,  ITOL,  IPASS  F0RMAT(4I5) 

IPRINT  -  PRINT  OPTION  FOR  DEBUGGING 

ICNVG  -  ITERATIVE  SOLUTION  TO  THE  CPHF  IS  CONVERGED  WHEN 
A*U  -  B  <  10**< -ICNVG).  DEFAULT=8. 

ITOL  -  IN  THE  ITERATIVE  SOLUTION  TO  THE  CPHF,  EXPANSION 
VECTORS  ARE  DROPPED  WHEN  THE  NORM  IS  <  10**(-IT0L). 
DEFAULTS  7. 

IPASS  -  NUMBER  OF  PASSES  ALLOWED  IN  THE  ITERATIVE  SOLUTION 
OF  THE  CPHF  EQUATIONS.  DEFAULT=2. 


-2-  IGOITR  F0RMAT(I5) 

=  0  USE  DIRECT  INVERSION  OF  A  IF  POSSIBLE, 

ELSE  USE  ITERATIVE  SOLUTION. 

<>  0  USE  ITERATIVE  SOLUTION  REGARDLESS  OF  SIZE. 
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SAMPLE  INPUTS 

Formaldehyde  SCF  optimization  in  an  ST0-3G  basis 


«  INPUT  «« 

FORMALDEHYDE  IN  ST0-3G  BASIS  AT  EXPERIMENTAL  GEOMETRY 


1000. 
CNV  2 

0  1 

0  10  10 

0 

OXYGEN 

GET  ST0-3G 

8. 

0.000000000000 

0.000000000000 

-2.286569522107 

CARBON 

GET  ST0-3G 

6. 

0.000000000000 

0.000000000000 

0.000000000000 

HYDROGEN 

GET  ST0-3G 

*  SCF  #### 

1. 

0.000000000000 

1.813400809181 

1.023872929647 

GROUND  STATE  SCF 


0 

10  0  0  0  0 

0  0  0  100 

5 

0 

0 

0 

1 

0 

2 

0 

0.0000000000 

1.0000000000 

«  TFOCK  ## 

GRSCF 

SCF  FIRST 

0 

«  NEVDER  # 

SCF 

FIRST  0 

XXXXX 

0  0 

I 
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A  subsequent  run  to  do  the  *A?  state  arising  from  a  b,,  -->  bj  excitation 
and  employing  the  old  vector,  would  have  the  same  input,  except  for  the  SCF 


«  SCF  #### 

EXCITED  SINGLET  STATE 
0  10  1  1  0  0 


000  100  00000000 


0.0000000000 

0.0000000000 

0.0000000000 

0.0000000000 


- 1 . 0000000000 
3.0000000000 
- 1 . oooooooooo 
1 . oooooooooo 


The  *2  state  would  have  the  following  input: 


tt  SCF  #### 

EXCITED  TRIPLET  STATE 
0  10  1  1  0  0 


000  100  00000000 


0.0000000000 

0.0000000000 

0.0000000000 

0.0000000000 


- 1 . oooooooooo 
-1 .oooooooooo 
-1.0000000000 
1.0000000000 


Fluorine  Azide  Cl  calculation  with  6-31G*  basis  set 


#  INPUT  ## 

FN3, Cl, 6-31G*, INTEGRAL  DATA 


1000. 

Cl 

0 

1 

0 

NITR0G1 

7. 

-1.0382783 

3.0050059 

1 

S 

6 

1 

4173.51 

0.00183477 

2 

627.458 

0.0139946 

3 

142.902 

0.0685866 

4 

40.2343 

0.232241 

5 

12.8202 

0.469070 

6 

4.39044 

0.360455 

2 

S 

3 

1 

11.6264 

-0.114961 

2 

2.71628 

-0.169117 

3 

0.772218 

1.14585 

3 

S 

1 

1 

0.212031 

1.0 

4 

P 

3 

1 

11.6264 

0.0675797 

2 

2.71628 

0.323907 

3 

0.772218 

0.740895 

5 

P 

1 

1 

0.212031 

1.0 

6 

D 

1 

1 

0.8 

1.0 

NITR0G2 

7. 

0.0 

1.2052621 

1 

S 

6 

1 

4173.51 

0.00183477 

2 

627.458 

0.0139946 

3 

142.902 

0.0685866 

4 

40.2343 

0.232241 

5 

12.8202 

0.469070 

6 

4.39044 

0.360455 

2 

S 

3 

1 

11.6264 

-0.114961 

2 

2.71628 

-0.169117 

3 

0.772218 

1.14585 

3 

s 

1 

1 

0.212031 

1.0 

4 

p 

3 

1 

11.6264 

0.0675797 

2 

2.71628 

0.323907 

3 

0.772218 

0.740895 

5 

p 

1 

1 

0.212031 

1.0 

6 

D 

1 

1 

0.8 

1.0 
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NITR0G3 

7. 

1.3921164 

1 

S 

6 

1 

4173.51 

0.00183477 

2 

627.458 

0.0139946 

3 

142.902 

0.0685866 

4 

40.2343 

0.232241 

5 

12.8202 

0.469070 

6 

4.39044 

0.360455 

2 

S 

3 

1 

11.6264 

-0.114961 

2 

2.71628 

-0.169117 

3 

0.772218 

1.14585 

3 

s 

1 

1 

0.212031 

1.0 

4 

p 

3 

1 

11.6264 

0.0675797 

2 

2.71628 

0.323907 

3 

0.772218 

0.740895 

5 

p 

1 

1 

0.212031 

1.0 

6 

D 

1 

1 

0.8 

1.0 

FLUORINE 

9. 

-0.2752066 

1 

s 

6 

1 

7001.71 

0.00181962 

2 

1051.37 

0.0139161 

3 

239.286 

0.0684053 

4 

67.3974 

0.233186 

5 

21.5200 

0.471267 

6 

7.40310 

0.356619 

2 

s 

3 

1 

20.8480 

-0.108507 

2 

4.80831 

-0.146452 

3 

1.34407 

1.12869 

3 

s 

1 

1 

0.358151 

1.0 

4 

p 

3 

1 

20.8480 

0.0716287 

2 

4.80831 

0.345912 

3 

1.34407 

0.722470 

5 

p 

1 

1 

0.358151 

1.0 

6 

D 

1 

1 

0.8 

1.0 

#  SCF  #### 
SCF  GROUND 

STATE 

FOR 

FN3 

00 

9 

0  0 

0 

0  000  0 

15 

0 

0.00 


-0.7114274 


-2.721322 


0  300  0  0  0  0  0  0 

1.000 


0.0 


0.0 


! 

J 

i 

i 

! 


i 

* 
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■7 

5 


#  DRT  #### 

DRT  GROUND  STATE  FOR  FN3 
0  2 
1  60 

15D0C1  45U0C1 


#  Cl  ####### 

0  40  0  1 


l 
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Formaldehyde  SCF  force  constant  calculation  in  double-zeta  basis 


#  INPUT  ## 

FORMALDEHYDE  IN  A  DZ  BASIS  AT  EXP.  GEOM. 

1000.  0  10  10 

CNV  2 

0  10 

OXYGEN  8.  0.000000000000  0.000000000000  -2.286569500000 

1  S  6 

1  . 781654000E+04  0. 203100000000E-02 

2  . 117582000E+04  0. 154360000000E-01 

3  . 273188000E+03  0. 737710000000E-01 

4  . 811696000E+(52  0. 247606000000E+00 

5  . 271836000E+02  0.611832000000E+00 

6  . 341360000E+01  0 . 241205000000E+00 

2  S  1 

1  . 953220000E+01  0. 100000000000E+01 

3  S  1 

1  . 939800000E+00  0. 100000000000E+01 

4  S  1 

1  . 284600000E+00  0. 100000000000E+01 
1  P  4 

1  . 351832000E+02  0. 195800000000E-01 

2  . 790400000E+01  0. 124189000000E+00 

3  . 230510000E+01  0. 394727000000E+00 

4  . 717100000E+00  0. 627375000000E+00 

2  P  1 

1  . 213700000E+00  0. 100000000000E+01 

CARBON  6.  0.000000000000  0.000000000000  0.000000000000 

1  S  6 

1  . 423260986E+04  0. 202900008298E-02 

2  . 634882019E+03  0. 155349997804E-01 

3  . 146097000E+03  0. 754109993577E-01 

4  .424973984E+02  0. 257120996714E+00 

5  . 141892004E+02  0. 596554994583E+00 

6  . 196659994E+01  0. 242516994476E+00 

2  S  1 

1  . 514769983E+01  0. 100000000000E+01 

3  S  1 

1  . 496199995E+00  0. 100000000000E+01 

4  S  1 

1  . 153300002E+00  0. 100000000000E+01 

1  P  4 

1  . 181557007E+02  0. 185339991003E-01 

2  . 398639989E+01  0. 115442000329E+00 

3  . 114289999E+01  0. 386206001043E+00 

4  . 359400004E+00  0.640088975430E+00 

2  P  1 

1  . 114600003E+00  0. 100000000000E+01 
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HYDROGEN  1.  0.000000000000  1.813400800000  1.023872900000 

1  S  3 

1  . 192406000E+02  0. 328280000000E-01 

2  . 289920000E+01  0.231208000000E400 

3  . 653400000E+00  0.817238000000E+00 

2  S  1 

1  . 177600000E+00  0. 100000000000E+01 


#  SCF  #### 

SCF  GROUND  STATE 

0  10  0000000  100  00000000 
5  0 

0  0 

1  0 

2  0 

0 . 0000000000  1 . 0000000000 

#  TFOCK  ## 

GRSCF  SCF  SECOND 
0 

#  NEVDER  # 

SCF  SECOND  0 
XXXXX  0  0 

#  DRT  ft### 

DRT  INPUT 

4  24 

5D0C1  9U0C1  1D0C3  3U0C3  2DOC4  4U0C4 
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4  **  * 
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